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In the open circular billiard particles are placed initially with a uniform distribution in their 
positions inside a planar circular vesicle. They all have velocities of the same magnitude, whose 
initial directions are also uniformly distributed. No particle-particle interactions are included, only 
specular elastic collisions of the particles with the wall of the vesicle. The particles may escape 
through an aperture with an angle 28. The collisions of the particles with the wall are characterized 
by the angular position and the angle of incidence. We study the evolution of the system considering 
the probability distributions of these variables at successive times n the particle reaches the border of 
the vesicle. These distributions are calculated analytically and measured in numerical simulations. 
For finite apertures 8 < n/2, a particular set of initial conditions exists for which the particles are 
in periodic orbits and never escape the vesicle. This set is of zero measure, but the selection of 
angular momenta close to these orbits is observed after some collisions, and thus the distributions of 
probability have a structure formed by peaks. We calculate the marginal distributions up to n = 4, 
but for 8 > 7r/2 a solution is found for arbitrary n. The escape probability as a function of n _1 
decays with an exponent 4 for 8 > n/2 and evidences for a power law decay are found for lower 
apertures as well. 

PACS numbers: 45.50.-j,05.20.-y,02.70.-c 



I. INTRODUCTION 

In the theory of dynamical systems, the study of the 
decay of simple hamiltonian systems is of much in- 

terest. These studies are an extension of earlier research 
concerning the closed version of such systems, where the 
main question is their ergodicity. The pioneering work 
in this area is centered on the system called Sinai bil- 
liard, a circular billiard with a smaller circular exclusion 
area in its interior Q, which was shown to be ergodic. 
Other two-dimensional ballistic billiards are known to be 
ergodic as well, such as the Bunimowitch stadium 0, 0] • 
One main interest in the studies of open billiards is the 
decay dynamics for long times. Bauer and Bertsch [lj 
found an exponential decay in an chaotic dynamics and 
a power law decay for a system with regular dynamics. 
These first results concerning the integrable case where 
later questioned by Legrand and Sornette pj, but it be- 
came clear that the difficulty in settling this question 
using numerical experiments is related to the high sensi- 
tivity of the results to initial conditions Q . In a more de- 
tailed simulational study for the chaotic two-dimensional 
Bunimovich stadium [4] , algebraic tails were found at suf- 
ficiently long times, but the weight of the algebraic tail 
tends to zero in the limit where the size of the aperture 
vanishes. 

In the classical circular billiard the (non-interacting) 
particles undergo elastic specular collisions with the wall. 
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Since two quantities are conserved (kinetic energy and 
angular momentum) , the system is integrable and there- 
fore non-ergodic [8| . A recent study of the open version of 
this system was undertaken by Vicentini and Kokshcncv 
|5| , and the algebraic long time of the survival probability 
was studied in detail, in the limit of a very small opening 
(weakly open billiard), using a coarse-grained description 
of the system. Here we study the problem making use of 
a probabilistic approach: starting from a random initial 
condition of the particles inside the vesicle, we calculate 
the joint probability distribution (in the angular position 
of point where the border is reached and the angle of in- 
cidence) at successive times the border of the vesicle is 
reached. This approach is not limited to small openings, 
actually it allows a full analytical treatment when the 
opening S is larger than n/2 (S — n/2 corresponds to a 
semi-circular vesicle). It will become clear below that for 
higher number of collisions with the wall, those particles 
which are close to periodic orbits which will never leave 
the vesicle will still survive, and thus a selection of these 
orbits will occur with increasing number of collision n. 
This selection of orbits according to the incidence angle 
(or, equivalently, to the angular momentum) was already 
found in the numerical experiments performed in Q , and 
emerges in a simple way in the probabilistic calculations 
below. The algebraic decay of the probability that the 
particle leaves the vesicle at the n'th collision is obtained 
analytically for large openings (5 > n/2). 

Let us define the problem. The initial position of a 
particle is distributed uniformly inside the circular vesicle 
of radius tq. The particle has a velocity of modulus vq, 
whose direction is also uniformly distributed. After some 
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FIG. 1: Movement of a particle, depicted at its initial position. 
This particle leaves the vesicle when it reaches the border for 
the second time. The indicated variables are described in the 
text. 



time, this particle reaches the border of the vesicle for 
the first time, at an angle 9\. If < S, the opening 
angle, the particle leaves the vesicle, otherwise it collides 
with the wall and performs the second straight part of 
its movement inside the billiard. If we use r$ as our unit 
of length and r Q / vq as our unit of time we may reduce 
the problem to the case of unitary radius and modulus 
of velocity. 

In Sec. |TT] we define the problem in more detail and 
present the analytical approach we adopted. The nu- 
merical study of the problem is described in Sec. IIII1 
and the simulational results are compared with the an- 
alytical ones and discussed. Further discussions and the 
conclusion may be found in Sec. IIVI 



II. 



DEFINITION OF THE PROBLEM AND 
ANALYTICAL APPROACH 



As discussed in Sec. HI we consider a particle inside a 
circular vesicle of unitary radius. The probability of the 
initial position, described by a radius r (0 < r < 1) and 
an angle 9q (— n < 9q < ix ), is uniformly distributed. The 
modulus of the initial velocity is unitary and the angle 
a (— 7r < a < 7r)between the position vector and the ve- 
locity is also uniformly distributed. Somewhat later, the 
particle reaches the border of the vesicle for the first time. 
If the angle 9\ where this occurs is such that — <5 < 9\ < S, 
the particle leaves the vesicle. Otherwise, the particle 
collides in a specular way with the border of the vesicle, 
continuing its movement until the border is reached for 
a second time at an angle 6>2- In figure Q] the trajectory 
of a particle is illustrated. 

The movement of the particle is deterministic, but the 
initial conditions are randomly distributed, and we are 
interested in the joint distribution of probabilities at the 



n'th time the border is reached. The initial joint distri- 
bution function is given by: 



p (r,a,9 Q ) = —2, 



(1) 



and this leads to the marginal distributions po(r) = 2r 
and po(a) — po{9) — 1/2-7T, which are consistent with the 
initial conditions stated above. Simple geometrical con- 
siderations allow us now to relate the variables 71 and t\ 
and 0\ at the first time the particle reaches the border 
with the original variables r, a, and 9 a . Notice that ti, 
the time of the first arrival at the border, is numerically 
equal to the distance between the initial position and the 
position of first arrival at the border. The angle of in- 
cidence 71 is defined in the interval — 7r/2 < 71 < 7r/2, 
and negative values correspond to particles moving clock- 
wise in the vesicle. With this convention, the signs of a 
and 7 are the same. The relations between the variables 
may be found applying trigonometric considerations to 
the initial segment of the trajectory illustrated in figure 
□ The result is: 
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The joint distribution at the first arrival at the border is 
related to the initial joint distribution as follows: 



Pi{~1i,h,Qi) = Po(r,a,9 ) 
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The jacobian may then be calculated using expressions 
H The result is: 
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and using the initial distribution [TJ we obtain the joint 
distribution at the first arrival at the border 
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From the joint distribution[S]the marginal distributions 
may then be calculated. The result for t\ is: 
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and the marginal distributions of the other two variables 
follow in a similar way: 
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As expected, the distribution of 9\ is uniform. The prob- 
ability that the particle leaves the vesicle at the first ar- 
rival at the border is given by: 



Pi 



Pi I 
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(9) 



Due to angular momentum conservation, the incidence 
angle will be the same at each subsequent time the par- 
ticle reaches the border. Also, the time interval between 
the first and n'th times at the border will simply be 
t' n = t n — t\ = 2(n — l)cos(7). We therefore consider 
this interval in the subsequent collisions and reduce our 
parameter space to the variables 7„ and 6 n . Integrating 
Pi(£ii7i>0i) over ti> we will get: 



pi {li, Ox) = ^cos 2 (7i). 



(10) 



Now we may relate the variables at the first and second 
times the border is reached: 



7i = 72, 

01 = 02 - 7T + 272. 
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The jacobian of this transformation is unitary, and thus 
the joint distribution at the second time the border is 
reached will be 

P2(72,0 2 ) = ^^/(01,*), (12) 

where the function f(6x,S) — 1 if |#i| > 5, and vanishes 
otherwise, so that only particles which did not leave the 
vesicle at the first time the border was reached do con- 
tribute. Actually 9\ in the right hand side of equation 
[T2] should be written as a function of 62 and 72 using the 
equations 1111 It is now easy to generalize this result for 
the n'th time the particle reaches the border: 
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We remark that in the absence of the aperture, the joint 
distribution will be the same at any collision, and there- 
fore for the closed circular billiard, the joint distribution 
[TOl and the marginal distributions [7] and [8] are the same 
for any value of n. For n > 1 and when the billiard is 
open, the joint distribution [T3] is not normalized, since 
particles leave the billiard through the aperture. 

Next, we may calculate the marginal distributions at 
the second time the particle reaches the border. The 
results are 
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and 
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One notices that the distribution in 9 2 is no longer uni- 
form if 8 > 0, displaying a maximum at 9 2 = 0. This may 
be explained qualitatively by the effect that due to the 
aperture, it will be more probable for a particle to have a 
positive than a negative horizontal component of the ve- 
locity after the first collision, since all particles which left 
the vesicle would have negative horizontal components if 
they had collided with the wall. Also, comparing the 
marginal distributions pi (71) and ^2(72), we notice that 
they differ just by a factor 1 — S/%, which accounts for 
the particles that left the vesicle at the first time the bor- 
der was reached. The distribution of t' 2 may be obtained 
from ,02(72): 



p 2 {t' 2 ) = 2p 2 ( l2 ) 
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where we defined t 2 = t' 2 /2, which is restricted to the 
interval [0,2]. The probability that the particle leaves 
the vesicle at the second time it reaches the border will 
be: 



P2 = 2 p 2 (0 2 )de 2 = 
Jo 

\[6(ir-6)+sm 2 {5)}. 



(17) 



As expected, this probability is larger than the value we 
would obtain for a uniform distribution in 9 2 , which is 
S(tt-S)/tt 2 . 

In general, to calculate the marginal distributions at 
the n'th time the particle reaches the border, we have 
to integrate the joint distribution 1131 For the incidence 
angle we obtain: 
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while the expression for the marginal distribution in 9 n 
is: 



1 [-tt/2 n-l 
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The distribution of the time t' n may be found generalizing 
equation [TBI The result is: 



Pn{t' n ) = 



T 2 I n (ln,S) 
(»-l)7T»0^f 



(20) 



where r„ = t' n /2{n — 1) and cos(7n) = t„. It may be 
useful to stress that all probability distributions are nor- 
malized with respect to the initial condition, and thus 
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FIG. 2: Graphical representation of IlTJi 1 f° r n = 4 

and 5 = 7r/10 as a function of #4 and 74. The product vanishes 
in the shaded areas and is equal to one in the complement. 
The horizontal segments at 71 — n/10 and 71 = correspond 
to initial conditions of a five-point star and diagonal peri- 
odic orbits, respectively, as discussed in the text. Particles 
with these initial conditions never leave the vesicle. The cir- 
cle indicates a particular crossing point of borders of stripes 
mentioned in the appendix. 



the area under each marginal distribution is an decreas- 
ing function of n, since particles are leaving the vesicle 
at each time the border is reached. 

To calculate of the integrals /„ and J n defined above, 
we have to identify the region in the (9 n ,2j n ) plane where 
at least one of the / functions in the product vanish. The 
region where a particular function f(9 n —j,8) vanishes is 
a set of j + 1 parallel stripes, as may be seen in figure 
[2J In this figure, which corresponds to n = 4, the two 
wider stripes correspond to (6*4, 274) values such that 
I #3 1 < 8, the three stripes of intermediate width corre- 
spond to 1 02 1 < 5 and points inside the four narrowest 
stripes are such that \6i\ < 6. The expressions for the 
stripes are obtained by a recursive application of equa- 
tions QT] which relate the variables at successive times the 
particle reaches the border. 

The integrals /„ and J n are continuous functions of 
their variables, but may display discontinuities in their 
derivatives. Their expressions are different in distinct 
ranges of their arguments. Although the calculation of 
these functions is rather straightforward, the identifica- 
tion of the ranges of integration becomes harder as the 
number of the collision n grows, the number of different 
expressions increases quite fast with n. In the general 
case, we performed the calculations up to n = 4. For 
8 > vr/2, however, the calculation of the marginal distri- 
butions may be accomplished for any value of n. Some 
details of these calculations and the results are shown in 
the appendix. 

Using the results for the marginal distribution p% we 
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The probability that the particle leaves the vesicle at its 
fourth arrival at the border is 
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Finally, using the marginal density for general n > 2 
obtained for 5 > n/2 in the appendix, we obtain: 
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III. DISCUSSION OF THE RESULTS AND 
COMPARISON WITH SIMULATIONS 

Initially, we show results for pi(ti) obtained with a 
simulation of 10 11 particles starting at the initial con- 
dition, that is, initial position inside the unitary circle 
with uniform probability density and the direction of the 
initial velocity also uniformly distributed. In figure [31 
the simulational results for the time of first arrival t\ 
are shown and compared with the analytical expression 
([6"]). In the same figure the numerical results for pi(ji) 
are compared with expression ([7]). It is apparent that 
the numerical results are compatible with the analytical 
curves. 

Although we realized many simulations starting from 
the initial condition, most of the numerical results shown 
here were obtained starting the simulation at the first 
collision of the particle with the wall (n = 1). With 
the same computational effort, these simulations led to 
much better results than the ones starting from the initial 
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FIG. 3: Results from simulations for pi(ti) (graph a) and 
Pi(ti-) (g ra P n b). The numerical results are indicated by cir- 
cles and the full lines are the corresponding analytical curves 
(equations [6] and [7] respectively) . 



condition, since only two variables define this condition 
now, instead of the three we have if the simulation starts 
at the initial condition. The numerical procedure followed 
in these simulations for a particular value of S was: 

1. Ng of values for 8±, uniformly spaced in the interval 
[— 7T, — S] U [S, 7r] are chosen. 

2. For each value of 8i, N 7 values for the incidence 
angle 71 are generated randomly with the proba- 
bility density [71 and the subsequent loci of arrival 
of the particle are calculated for n = 2, 3, . . ., up to 
a maximum number n max . 

3. At the n'th arrival at the border, the values of 9 n , 
-f n , and t' n are recorded, so that estimates for the 
marginal probability densities of these variables are 
obtained. If —5 < 9 n < 5, the counter of the num- 
ber of particles leaving the vesicle at the n'th arrival 
at the border is increased by one and the simula- 
tion of a new particle starts. If n — n max , the 
simulation of a new particle is started, otherwise 
we proceed to step n + 1. 

The results shown below were obtained with the choices 
N e = 10 5 , A 7 = 10 5 , so that a total of 10 10 particles are 
considered in each simulation. We chose n m ax = 30. 

In figure 0] we show simulational and analytical results 
for the marginal distribution p(j) when 6 — 0.5. One no- 
tices that the agreement between simulation and theory is 
very good for the first four values of n, much better than 
the one observed in figure [3] This is due to the reduc- 
tion of the number of initial conditions mentioned above. 
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FIG. 4: Results for pni'Jn) for n = 2, 3, ... 8. For the same 
value of 7, p„ is a non-increasing function of n. For n = 
2, 3, 4 the numerical results are indicated by circles and the 
full lines are the corresponding analytical curves (equation 
1181) . For larger values of n, the lines correspond to results of 
simulations. Values of n increase downwards. 



Another feature which is apparent in these data is that 
as n grows, peaks appear in the distribution at specific 
values of 7. The first peak to appear is located at 7 = 
and is already apparent in the analytical expression for 
Piijij- These peaks are associated to periodic orbits in 
the vesicle, already discussed in [8| . Between two arrivals 
at the border, the angle 9 changes as A8 — tt — 2-f. If the 
orbit closes after p collisions and the particle performs q 
turns in the vesicle, we have pA9 = 2nq. Thus, the angle 
7 Pi9 for this periodic orbit is: 

P- 2 g , 07 x 

where p and q have no common divisors, p > 2q and we 
considered only particles moving counterclockwise (A0 > 
0). The periodic orbits may therefore be associated to the 
rational numbers q/p. For example, p= 2,q= 1, corre- 
sponds to 7 = 0: the particle moves along a diameter. A 
polygonal orbit with p vertexes corresponds to (p, 1). A 
simple star with 5 vertices's is described by (5, 2). 

In an open vesicle, only those periodic orbits with 
7r/p > 5 will last, and the peak will develop when n > p. 
This may be simply understood if we consider that a par- 
ticle in a periodic orbit characterized by an index p will 
reach the border of the vesicle at p equally spaced points. 
If the angle between two fist neighbors is larger that 25, 
for appropriate values of Q\ the particle never reaches the 
aperture and therefore remains inside the vesicle forever. 
Of course these orbits are a set of zero measure in the 
initial conditions. The four peaks visible in figure [H at 
non-vanishing 7 may be identified as (5,2) (7 = tt/10), 
(3, 1) (7 = tt/6), (4, 1) (7 - tt/4), and (5, 1) (7 = 3tt/10). 
The hexagon (6, 1) still fulfils the condition above for the 
aperture S = 0.5 (n/6 > 0.5), but since it is quite close 
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FIG. 5: Results for p n (7n) for n — 4, 5, 6 and S — 1.6. For 
n = 6 the numerical results are indicated by circles. Full lines 
are the analytical curves (equation [T8]) . 




FIG. 6: Results for p n (0 n ) for n = 2, 3, 4 ... 8 and 8 = 0.5. For 
n — 2,3,4 the numerical results are indicated by circles and 
the full lines are the analytical curves (equation [T^J). For n > 
4, only simulational results are shown. Values of n increase 
downwards. 



to the limit the corresponding peak at 76,1 = 7r/3 is very 
small and hard to identify in the plots above. It is visi- 
ble, however, in the results for larger values of n. Once 
a peak appears, it will narrow with increasing n. This 
may be explained noting that, of the particles with the 
corresponding angle 7 Pjl? , some will leave the vesicle, but 
only when n < p. Thus, the maximum value at the peak 
is stationary for n > p. However, particles with a slightly 
different incidence angle still leave the vesicle for n > p. 
As n — > 00 the peak shrinks to just one point. The condi- 
tion on the initial value of 9 which assures that a particle 
with 7 = 7 Pj9 will never leave the vesicle is £\9±\ > 8, 
for I = 1, 2, . . .p and reducing each angle to the interval 
[— 7r, 7r[. In summary, the periodic orbits correspond to 
horizontal segments in the (#1,71) space of initial condi- 
tions. The initial conditions for a periodic orbit labeled 
as (q,p) in a vesicle with aperture 8 are p segments at 
71 = 7 P .g, with length A0i = 2(it/p — 8) and endpoints 
\9\ at (S + 2en/p,8 + £iT/p + 89), for I = 0, 1, 2, . . . ,p/2. 
In figure [5] these segments were drawn for the star orbit 
(5, 2) and for the diagonal orbit (2, 1). 

For 8 > 7r/2, all periodic orbits are suppressed, and 
we are able to calculate the marginal distribution for any 
n. In figure [5] results for 8 = 1.6 are shown for n — 
4, 5, 6. The absence of peaks related to any periodic orbit 
is apparent, and for n — 6 the agreement between the 
analytical result (line) and the simulations (circles) may 
be appreciated. In this case, all particles with \j\ < ((n+ 
1)8 — 2tt)/(ti — 1) leave the vesicle at the n'th time the 
border is reached or before. 

Marginal distributions for 9 are shown in figure [5] and 
compared with the analytical results for n < 4, also for 
the case 8 — 0.5. Again some peaks are observed, but 
they are much broader than the ones in p("f). If we con- 
sider only a periodic orbit with p vertexes, we would 
have p equally spaced rectangular peaks, of width A9i 
defined above. Thus, as p grows approaching the limit- 



ing value n/8 the peaks become narrower. Therefore, for 
8 = 0.5 the narrowest peaks belong to hexagons, and it is 
this peak we may identify for largest values of 9, located 
close to 9 = 2.61. For lower values of 9 we have the su- 
perposition of several peaks and identification is harder. 
However, for larger n a step like structure develops at 
smaller 9, as expected. Although we do not show results 
here for 8 > tt/2, where we obtained p n (9 n ) analytically 
for any n, as expected smooth curves are found, since no 
periodic orbits survive in this case. All particles with 
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leave the vesicle at the n'th arrival at the border or be- 
fore. 

In figure [7] we show results for the marginal distribu- 
tions in the time t' n between the n'th and first arrival 
of the border. Again a good agreement between analyt- 
ical end numerical results is found for n = 2,3,4. As 
expected, peaks develop at values which correspond to 
surviving periodic orbits, as may already be seen in the 
simulational results for larger values of n. The largest 
peak located at t„ = t' n /2(n — 1) = 1 corresponds to the 
diagonal periodic orbit (2,1), the other periodic orbits 
originate peaks at lower times. For 8 > n/2 again results 
for general n may be found and their agreement with 
simulational results is very good. Since just polygonal 
paths 9=1 survive after some collisions in this case, the 
maximum time a particle may spend inside the vesicle is 
2(ir — 8). This result may be obtained from the general 
expression of p n (t' n ) in the limit n — > 00. 

Figure [8] shows results for the probability p n that a 
particle leaves the vesicle as the border is reached the 
n'th time, as a function of the aperture 8 for the cases 
where analytic results are available. A good agreement is 



FIG. 7: Marginal density of probability of the time interval 
since the first arrival of the border, as a function of r n , for 
n = 2,3, ... ,8 and S — 0.5. Full lines for n < 5 correspond 
to analytical results (equation I20[) , and the circles are results 
from the simulations. For higher values of n only numerical 
results are presented. Values of n increase downwards. 




FIG. 8: Plots of p„ as a function of S for n = 1, 2,3,4. Re- 
sults from simulations are represented by circles and the lines 
are analytic curves (equations l9l 117112111221 and 1231251 respec- 
tively) . 



FIG. 9: Plots of p n as a function of n for S — 0.1 (circles), 
S = 0.5 (squares), and 6 = 1.0 (diamonds). The thick lines 
are exponential functions discussed in the text. 



circular billiard 5 = the distribution p{9) is uniform, 
so that, at least for small n, as S increases the range 
of values of the probability density becomes larger. As 
mentioned above, P2 is always larger than the value for 
an uniform distribution, but in general we observe that 
for an intermediate range of n, the exponential values 
are larger than the numerical results, but at sufficiently 
large n they again become smaller. For 5 > 7r/2, we may 
expand the expression 1261 in l/n and find the asymptotic 
power law behavior: 
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so that in this case we may assure that p n is larger than 
Pn jU for sufficiently large values of n. It remains an open 
question if this conclusion extends to S < 7r/2, but our 
numerical results are consistent with this possibility. For 
n = 3,4, p n is larger than p n<u for S in the range ]0,po[ 
and smaller in the range ]po, 7r[, where pa = tt/2 for n = 3 
and po = 0.919725 . . . for n = 4. 



IV. CONCLUSION 



found between numerical and analytical results. In figure 
|5]we present numerical results for function of n for 

some values of 5. We also included exponential functions 
which would be obtained for the escape probability if the 
distribution p{6) would be always uniform, which are: 

Pn,u=(l--) -■ (28) 
\ IT J IT 

First, we notice that the difference between the expo- 
nential decay and the results is lower for small values of 
8. This may be explained recalling that for the closed 



In our study of the open circular billiard, it is natural 
to obtain the distribution of the quantities of interest as a 
function of the collision number n. In other approaches, 
as usual, the time appears as the variable in terms of 
which the probability distributions are calculated. Al- 
though we did discuss somewhat the distribution of times 
at the first collisions, we did not enter in much detail in 
this question here. This is the subject of ongoing work. 
In the simulations, we also did not record results beyond 
30 collisions with the wall, so that the question of the 
asymptotic decay of p n for S < tt/2 was not discussed in 
detail. As stated above, we have indications of a power 



law behavior, but the exponent may be different from 
the value 4 found exactly for larger apertures. One may 
imagine that the existence of orbits in which the parti- 
cles are trapped forever inside the vesicle for 6 < tt/2 
may lead to a slower decay, with exponents smaller than 
4 in this case. It is true, as may be seen above, that these 
orbits are a set of zero measure in the initial conditions, 
but particles in a finite region close to those orbits may 
take a long time to exit. 

As was shown, the calculation of marginal distributions 
in the general case becomes more difficult for increasing 
n and finite S < tt/2. There is, however, a simplification 
when 8 <C 1, that is, if we consider the first terms in the 
expansion of the distributions in powers of 8. This is the 
range where other approaches, such as the one in are 
effective. 

There is a promising possibility to access the asymp- 
totic behavior at long times if we recall that, after a cer- 
tain number of collisions with the wall, the particles re- 
maining in the vesicle will have incidence angles close to 
the values defined for the surviving periodical orbits. As 
another illustration of this, in figure [TUl numerical results 
for the distribution of the time interval t' 2 Q = t 2 Q — ti 
are depicted in a logarithmic scale for a simulation with 
8 = 0.5. The most important peak at t' 20 = 38 corre- 
sponds to the vicinity of the diagonal orbit (2, 1). The 
other peaks in order of decreasing values of t' 20 are the 
same mentioned in the discussion of figure [4] above, but 
now the small peak corresponding to the hexagon (6, 1) is 
visible centered at t' 20 — 38cos(76,i) = 38cos(7r/3) = 19. 
The small continuous distribution at low times corre- 
sponds to particles in high p orbits (whispering gallery 
modes), which still did not reach the aperture after 19 
collisions. As discussed above, this distribution is located 
at times smaller than 2(7r — 5) and does not contribute at 
larger times. We are presently working on an approach 
where the neighborhood of the periodic orbits is treated 
in a approximate way, which is should be appropriate for 
the discussion of the long time limit. 
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APPENDIX A: CALCULATION OF /„ AND J„ 

Both I n (j n , 8) and J n (9 n , 8) are even functions of their 
first argument, so we restrict the calculations below to 
non-negative values of these variables. As may be seen in 
figure [U the domains of integration are split into regions 
where the product of / functions is equal to unity. This 
regions change as the variables of the marginal distri- 
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FIG. 10: Probability density of the time interval t'20 obtained 
from simulations with 8 — 0.5. The peaks which correspond to 
the neighborhood of surviving periodic orbits identified in the 
text are apparent. The broad peak at low times corresponds 
to whispering gallery modes still present in the vesicle. 




FIG. 11: Regions of the domain of the arguments where 
^3(7n, 8) is given by different expressions. The function I3 
vanishes in the triangular region labeled by 0. 



butions are varied, so that these distributions are given 
by different expression in different areas of the domain 
of their arguments. It is convenient to use the variable 
<fi n = 2"f n in the functions I n . For I3 these areas are 
shown in figure [TT1 

For a particular value of the argument 03, I3 is equal 
to to the sum of the segments in white regions of the 
corresponding horizontal line in figure if we erase the 
narrowest stripes which correspond to n = 4. For 8 < 
tt/2, this total length is equal to 2-7T — 48 if 03 < ir — 28 
and to 7T + 03 — 28 otherwise. In general, the expressions 
for 73 in the triangular regions depicted in figure [3] are 

/3 1 = 2n-4S, (Al) 
if = 7T + 3 -2(5. (A2) 



jc/2 



jc/2 



FIG. 12: Regions of the domain of the arguments where 
Js(7n,5) is given by different expressions. The equations of 
the limiting lines are 83 = 6, 83 = 35, 83 = 2tt — 35, and 
83 = 5 — 2-7T. J3 vanishes in the region labeled by 0. 



FIG. 13: Regions of the domain of the arguments where 
5) is given by different expressions. The equations of 
the limiting lines are 4>4 — 5 and ^4 = 7r — 25. I4 vanishes in 
the region labeled by 0. 



In these calculations, the crossings of the borders of the 
stripes determine the limits of the region in the argu- 
ments where each expression is valid. These borders have 
the general form 



(n - 2£)ir ± S, 



(A3) 



where the plus sign is for the upper border and the minus 
sign for the lower border of the stripe and < I < n 
labels the stripes in each order n. As an example of such 
a crossing of borders, we have the point indicated by a 
circle in figure [2 which corresponds to the crossing of 
the lines (n = 3,£ — 0, — ) and (n — 2,1 = 0, +) and is 
located at 9 — 3<5, <fi = it — 25. 

For the calculation of J„ , each vertical segment in the 



white region of the diagram[211imited by 0_ and 
tributes with g(4>+) — .9(0-) to the function J n 
integrating equation (flUj) we get: 
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(A6) 
(A7) 
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(A9) 



For n = 4, the function 74(74, 5) is given by three dif- 
ferent expressions in the triangular regions indicated in 
figure 1131 vanishing in the fourth region in the space of 
its arguments. The expressions are: 



9{4>) = -[0 + sim»]. 



(A4) 



H 
It 



2(tt - 04 - 25) 
2(tt-3J) 
2(04 -<5). 



(A10) 
(All) 
(A12) 



We may then obtain the expressions for J 3 . The diagram 
indicating the region where each expression is valid is 
shown in the figure [T2l The expressions are: 



Finally, in figure[T3]the domain of the arguments of the 
function J±(6±,S) is depicted, with a total of 18 regions 
where different expressions are valid. The expressions 
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FIG. 14: Regions of the domain of the arguments where 
J4,(di,8) is given by different expressions. The equations of 
the limiting lines are O4 — 5, 64 — 35, 64 — 55, 64 — 2n — 55, 
64 = 7T - 25, 64 = tv - 5, 64 = 2tt - 35, 9 4 = -tt + 25, and 
84 = 55 — 2ir. J 4 vanishes in the region labeled by 0. 
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The extension of these calculations to larger values of n 
becomes increasingly difficult and will not be presented 
here. However, for n > 2 and <5 > 7r/2 a major sim- 
plification occurs and it is easy to perform the calcu- 
lations of and J n +ii since in the diagram corre- 
sponding to figure [5] in this case the white area is re- 
duced to only two triangles and two quadrangles. The 
vertexes of the first triangle are located at 9 = (f> = — 7r, 
= —it, (f> = — 7r(l — 1/n), and 9 = —5, <fr = — tt. The sec- 
ond triangle may be obtained by transforming 9^—9 
and (/>—>—(/). The vertexes of the quadrangle are lo- 
cated at 9 = n,4> = — 7T, 9 = n, (j) = — 7r(l — 1/n), 
6 = 6, <j) = -7T, and 9 = (5{n + 1) - 2n)/(n - l),<f> = 
— ((n — 3)7r + 2S)/(n — 1). The second quadrangle may 
be obtained by the same projection operation described 



above. We then obtain 

I n+1 = O,forO<0 w+1 < fr-ff*+ a * ; 
Z„ +1 = (3-n)7T-2<5+(n-l)^„ + i, 



otherwise. 



(A32) 



For the function related to the marginal distribution 
in we get the result J„+i = for < 6* n+ i < 

((n + l)<5-27r)/(n- 1) and: 
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, for S < 9 n+1 < tt. (A34) 
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